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A new search for features in the primordial power 
spectrum 
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ABSTRACT 

We develop a new approach toward a high resolution non-parametric reconstruction of 
the primordial power spectrum using WMAP cosmic microwave background temper- 
ature anisotropics that we confront with SDSS large-scale structure data in the range 
k ~ 0.01 — O.lftMpc -1 . We utilise the standard ACDM cosmological model but we 
allow the baryon fraction to vary. In particular, for the concordance baryon fraction, 
we compare indications of a possible feature at k ~ 0.05 /iMpc -1 in WMAP data with 
suggestions of similar features in large scale structure surveys. 

Key words: cosmology: early universe - cosmic microwave background - large-scale 
structure of Universe - cosmological parameters — methods: data analysis. 



1 INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe (WMAP) ex- 
periment has recently provided a high resolution full sky 
cosmic microwave background (CMB) map that provides 
an unprecedented opportunity to probe the structure and 
the contents of the universe (Spergel et al. 2003). The cor- 
responding angular power spectrum has been used to set 
constraints on cosmological parameters. In the framework of 
ACDM cosmologies, such constraints are strong if the initial 
conditions are parametrized by primordial power spectra in 
the form of power laws with or without running spectral in- 
dex. However parameter estimation can be weakened and bi- 
as ed if more complex fo rms of initial conditions are adopted, 
cf. lBucher et alJ 12004ft and lBlanchard et alJ fcOOSt . We ex- 
plore another such possibility here. 

Recent papers from the Sloan Digital Sky Survey 
(SDSS) team report measurements of the matter power spec- 
trum and the cosmological parameters that give a best fit to 
its shape (Tegmark et al. 2003a,b; Pope et al. 2004). In par- 
ticular Pope et al. (2004) finds that for the for a value of the 
matter density, £l m h = 0.264 ± 0.043 and baryon fraction, 
fb = Qb/n m = 0.286 ± 0.065 at 1 a. These determinations, 
in particular that of the baryon density, are in marginal con- 
flict with the data from both WMAP (Spergel et al. 2003) 
and primordial nucleosynthesis (BBN) (Cyburt et al. 2003; 
Cuoco et al. 2003), which, together with the Hubble con- 
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stant (Freedman et al. 2001) determination, prefer a lower 
baryon fraction, ~ 0.2 when a power law primordial power 
spectrum in a ACDM cosmology is adopted. 

Such a discrepancy might be due to the presence of in- 
trinsic features in the power spectrum that affect the SDSS 
(and other survey) data, such as for example a dip and a 
bump at k ~ 0.035 /iMpc -1 and k ~ 0.05 /iMpc -1 (see also 
Atrio-Barandela et al. 2001; Barriga et al. 2001). It is useful 
to explore the implications of such a possibility, even though 
the data at present is not compelling, because very similar- 
features in the matter power spectrum have been detected in 
independent large-scale structure surveys, notably the 2dF, 
Abell/ACO cluster and PSCz surveys (Einasto et al. 1997; 
Percival et al. 2001; Miller & Batusky 2001; Hamilton et al. 
2000). In this work, we show that, if the power law hypoth- 
esis for the shape of the initial power spectrum is relaxed, 
a still good or even better agreement between the different 
data sets can be obtained. Our strategy is principally com- 
prised of two steps: a high definition reconstruction of the 
primordial power spectrum from WMAP data (for given sets 
of cosmological parameters) and a convolution of the initial 
power spectrum with the SDSS window functions to deter- 
mine the statistical agreement between the two data sets. We 
choose to study only these two data sets in a limited range 
of k— space in order to make a preliminary exploration of 
the possible influence of primordial features. We find that 
for reasonable choices of the cosmological parameters, the 
WMAP-reconstructed matter power spectrum shows oscil- 
lations that are very similar to those of SDSS once the convo- 
lution is performed. This results in an improved concordance 
for the inferred baryon fraction and at the same time is sug- 
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gestive of a deviation from a power-law primordial power 
spectrum. 



2 METHOD 

The temperature CMB angular power spectrum is related to 
the primordial power spectrum, P°, by a convolution with a 
window function that depends on the cosmological parame- 
ters: 



Ce = 4tt 



A 2 e (k)P°(k)f 



WP U 



(1) 



where the last term is in matrix notation and represents a 
numerical approximation to the integral, calculated using a 
modified version of CMBFAST (Seljak & Zaldarriaga 1996). 
The Ce spectrum then contains information about both the 
cosmological parameters and the initial power spectrum. 

Various papers have presented reconstructions of the 
power spectrum from WMAP data (e.g. Bridle et al. 2003; 
Mukherjee & Wang 2003; Matsumiya et al. 2003; Shafieloo & 
Souradeep 2004; Kogo et al. 2004; Hannestad 2004). Most 
of these have confronted an a priori parametrised power 
spectrum with the data to obtain inform ation about its 
shape and amplitude. Following the spirit of lGawiser fc Silkl 
( 1998) and lTegmark fc Zaldarriagal (120021) . our power spec- 
trum is not described by a small set of parameters that 
incorporate features, but is finely discretized in fc-space 
(see also Shafieloo & Souradeep 2004; Kogo et al. 2004). 
Our work represents the first attempt to reconstruct the 
power spectrum at high resolution in the full range k ~ 
0.01 — 0.1 /iMpc -1 together with an estimation of the error 
covariance matrix. Our findings are cons istent with t he re - 
sults obtained with a different method bv lKoeo et al] (120041) 
that are limited to k < 0.045 /iMpc -1 , corresponding to 
I < 430. 

Generally speaking, it is not possible to invert the ma- 
trix W to solve for the power spectrum, because each Ce 
embraces information about a limited range in fc-space. How- 
ever the inversion is feasible under some assumptions. Basi- 
cally the lack of information can be remedied by the intro- 
duction of priors, such as for example the requirement of a 
certain degree of smoothness in the solution. To be specific, 
we consider a solution of the form 



P = MCT MA ^ 

with an error covariance matrix given by 
E = [W'N^W + eL t L]~ 1 . 



(2) 



(3) 



WMAP 

e are 



Here: M = [W'A -1 W + elSL] 1 W t N~ 1 ; C, 
the mean values of WMAP data; N is the WMAP co- 
variance matrix (Verde et al. 2003) and L is a discrete 
approximation to the first derivative operator. Our ap- 
proach follows the spirit of a similar solution proposed in 
iTegmark fc Zaldarriagal (120021) : however we replace the iden- 
tity operator by L. The parameter e regulates the degree of 
solution smoothness, since the derivative operator L acts in 
such a way to minimize the norm of the solution derivative. 
This solution can be effectively thought of as a linear least 
squares solution modified to accommodate smoothing. 



In order to select a value for the parameter e, we calcu- 
late the angular power spectrum Ce resulting from P. Then 
we form the expression 



X = {Ce 



■Ce WMAP ) t N~ 1 (C l -C e WMAP ). 



(4) 

It is reasonable to fix the parameter e to render the \ 2 m 
the previous equation equal to the number of WMAP tem- 
perature data points, namely 899. In other words, a solu- 
tion is chosen that gives an acceptably good fit to the data. 
Moreover, as will be explained below, this represents a con- 
servative choice since the smoothing decreases the effective 
number of degrees of freedom. 

To test the method, we reconstruct the power spec- 
trum by using angular power spectra calculated from var- 
ious power spectrum shapes with and without added er- 
rors. We obtain very good results, except when there are 
features of extension comparable to our fc-space gridding 
Afc, which we took to be approximately 10 -4 in the range 
0.01 — 0.1 Mpc -1 . The integration range is typically between 
10 -5 and 0.3 Mpc -1 and comprises about 1500 points. We 
also verify that the signals detected in the actual WMAP 
data (see below) do not depend on discretization issues by 
testing different resolutions. Our method yields a fc-space 
resolution that is limited only by computing resources, at the 
price however of introducing correlations between neighbor- 
ing points that we take fully into account into our statistical 
analysis. 

In analogy with ordinary regression it is possible to de- 
fine the effective number of degrees of freedom caused by 
smoothing by the expression Trace (WM). It is easy to rec- 
ognize that without smoothing this quantity would have 
been equal to the number of bins in fc-space. We have found 
that our choice of the parameter e implies the effective us- 
age of approximately 45 constraining informations out of the 
more than 1000 amplitudes that specify the power spectrum. 

The procedure just described must be repeated for each 
given set of cosmological parameters. We proceed by first 
evolving the primordial power spectrum to redshift zero by 
multiplying it by the appropriate transfer function calcu- 
lated by CMBFAST. Then we convolve the derived matter 
power spectrum with the SDSS window functions given by 



Tegmark et al . (2003a). Finally, we evaluate 
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where p SDSS an d N SDSS represent the SDSS mean data 
and (diagonal) covariance matrix (Tegmark et al. 2003a), P 
and E are the SDSS-convolved WMAP-reconstructed mat- 
ter power spectrum and covariance matrix (calculated by 
error propagation) and b is the bias parameter. We follow 
the common practice of considering a constant bias which 
we set at the optimal value. 



3 RESULTS 

In performing the \ 2 analysis indicated by Eq. (J5J, we re- 
strict the computation to the SDSS experimental data cor- 
responding to the 11 points with lower fc-values, such that 
k < 0.1/iMpc _1 , in order to avoid effects caused by the 
limitation in fc-range that arises since we are extracting in- 
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Table 1. Models, inferred optimal bias, and \ 2 comparison be- 
tween WMAP constrained matter power spectra and the SDSS 
data for a selection of the first 11 points with lowest wavenumber, 
corresponding to 10 degrees of freedom. 



Model 


fb 






h 


b 


X 2 /dof 


1 


0.171 


0.194 


0.024 


0.72 


1.04 


4.7/10 


2 


0.286 


0.264 


0.054 


0.72 


0.97 


26.4/10 


3 


0.155 


0.215 


0.023 


0.695 


1.10 


3.3/10 


4 


0.200 


0.209 


0.030 


0.72 


1.11 


3.6/10 


5 


0.143 


0.194 


0.020 


0.72 


1.00 


6.4/10 



formation from the CMB power spectrum and the quality 
of the WMAP data degrades significantly above £ ~ 800. 

We select a few choices of the cosmological parame- 
ters, for which we perform the reconstruction. The pres- 
ence of features in the power spectrum is associated with 
the deviati ons from a power law angular power spectrum 
noticed by ISpergel et alJ i2003T) . These authors mention 
that the first release of data did not include effects con- 
tributing roughly 0.5-1% to the power spectrum covariance, 
mainly due to gravitational lensing, beam uncertainties and 
non-Gaussianity in the noise maps. Therefore, to verify if 
they are real, a cross-check with an independent dataset is 
needed. The SDSS data serve such a purpose. 

In Table we list the models used, defined in terms 
of the baryonic and total matter content, Hubble constant 
h (in units of 100 km sec -1 Mpc -1 ) and bias b. The corre- 
sponding x 2 f° r a selection the first 11 SDSS data points 
with lowest wavenumber are also given and will be used for 
the discussion below. The optical depth r is set equal to the 
best WMAP case, 0.166. We have chosen to concentrate on 
a few cases that differ mainly with respect to the baryonic 
fraction, since this is the most interesting quantity to probe 
in light of the SDSS analysis that uses a power-law power 
spectrum. 

We label the models as: model 1, the WMAP mean 
best case (Spergel et al. 2003); model 2, the SDSS best case 
(Pope et al. 2004) (in marginal conflict with the previous 
one); model 3, close to the joint WMAP-SDSS best case 
(Tegmark et al. 2003b); model 4 and model 5, characterized 
by a relatively high and low baryonic content. The "best 
cases" we refer to were all obtained for a power-law power 
spectrum. The latest constraints from primordial nucleosyn- 
thesis give, at 1 a: uj b = h 2 fl b = 0.022 ± 0.002 (Cyburt et 
al. 2003; Cuoco et al. 2003), which excludes model 2 with 
cu b = 0.05 . 

Our first result is a reconstruction of the initial power 
spectrum shape from the CMB angular power spectrum with 
the associated error covariance matrix. Fig.Qshows the cor- 
responding initial spectrum reconstructed from the WMAP 
data and assuming the cosmological parameters of model 3, 
in the range 0.01 < k < 0.075 Mpc -1 . The effective range 
probed in t is approximately 140-800. By reconvolving the 
power spectrum back to multipole space, it is possible to 
note in Fig. |2] the action of smoothing and appreciate the 
accuracy of the inversion. Even if not far from a power law, 
small features are present in the initial power spectrum at 
all scales; we shall devote our attention mostly to the suc- 
cession of bumps and dips around 0.035 Mpc -1 . 

It is useful to implement Monte Carlo simulations to 
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Figure 1. Reconstructed initial power spectrum in units of 
2.95 X 10 -9 from WMAP Cf's with model 3 cosmological pa- 
rameters. The irregular middle curve represent the mean and the 
surrounding curves are obtained by summing and subtracting the 
square root of the diagonal elements of the covariance matrix. 
The smooth line is a reference power law power spectrum given 
by fc" s_1 , where n s ~ 0.99. 
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Figure 2. The noisy thin curve connects the data from WMAP. 
The reconstructed power spectrum convolved back to multipole 
space is given by the thick smooth line, that shows clearly the 
effect of smoothing. 



look for deviations and systematic biases in the reconstruc- 
tion. We generated 5000 random realizations around the an- 
gular power spectrum calculated from a power law primor- 
dial power spectrum assuming a gaussian distribution, which 
is reasonable in the probed multipole range. The case for 
model 3 is presented in Fig. |3] The averaged reconstruction 
doesn't show strong bias and deviations bigger that the one 
sigma band, obtained from the standard deviation respect to 
the averaged reconstruction, are also evident. Again, there 
are indications of possible features at k ~ 0.035 Mpc -1 and 
other wave numbers. 

Models 1 and 2 reconstructed matter power spectra are 
plotted in Fig. 21 Model 2, constructed to fit the SDSS data, 
is marginally excluded by WMAP data under the initial 
power law assumption: it requires too many baryons in order 
to fit the features (oscillations) of the SDSS reconstructed 
matter power spectrum. If we relax this assumption by us- 
ing our method and keeping the same baryon fraction, we 
see that the WMAP Gt features translated into k-space do 
not match the SDSS features (x 2 = 26/10). Baryonic oscil- 
lations and CMB features do not converge to r eproduce the 
SDSS behaviour. As noted in lPope et alJ <l2004t) . oscillations 
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0.01 0.015 0.02 0.03 0.05 0.07 

k [ Mpc -1 ] 

Figure 3. Results from 5000 Monte Carlo realizations drawn 
form a power law power spectrum, given by fc" a — 1 where n s = 
0.99, with model 3 cosmological parameters. The thin smooth 
line in the center is the exact initial power law and the dashed 
curve is the averaged reconstructed initial power spectrum. The 
two symmetric lines represent the standard deviation from the 
averaged reconstruction. No strong systematic bias appears in 
the reconstruction and possible candidates for features are given 
by the outliers respect to the one sigma band. The plot is in units 
of 2.95 X 10" 9 . 
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induced by the prefered baryon fraction value of CMB are 
out of phase with SDSS "oscillations". 

Conversely, assuming a power law initial power spec- 
trum, model 1 gives a poor fit to CMB data (Spergel et al. 
2003) and a good fit to the selected SDSS data (\ 2 =6/10). 
When we allow for a free shape for the initial power spec- 
trum, using our method, model 1 becomes by definition a 
good fit to WMAP Ci's and decently fits the 11 points se- 
lected from SDSS (x 2 = 5/10). By taking into account the 
features in the CMB we improve not only the fit to the CMB 
but also to the SDSS data. 

However we found other models that can achieve a bet- 
ter fit. If we consider the cosmological parameter of model 
3 (close to the best WMAP+SDSS fit), our reconstructed 
matter power spectrum gives a better fit (\ 2 ~ 3.3/10, see 
Fig. In doing this, we show that the same model could 
fit both the WMAP GYs and SDSS matter power spectrum 
shape. This last point is our principal result, and suggests 
that an acceptable compromise that simultaneously fits both 
WMAP and SDSS may require a change respect to the con- 
cordance model, such as deviations from a power-law initial 
power spectrum. 

We then investigated the effect of varying the baryonic 
fraction in model 3 and find that features from the CMB 
angular spectrum must be combined with the appropriate 
amount of baryons in order to fit the SDSS "oscillations" 
(see Fig.QjJl. Model 5, with a low baryonic content, is a rela- 
tively bad fit because the WMAP features seem almost out 
of phase with the SDSS features. While model 4 leads to 
a better agreement since the shapes of the SDSS k ~ 0.035 
and 0.05 foMpc -1 features are mimicked reasonably well. We 
furthermore note that model 3 falls in the observational 
range of both h (0.72 ± 0.07, Freedman et al. 2003) and 
n b h 2 (0.022 ±0.002, Cyburt et al. 2003; Cuoco et al. 2003). 

Finally, looking at Fig. [(J suggests that the feature seen 
at k ~ 0.05 /iMpc" 1 in the SDSS is well reproduced by the 
reconstructed WMAP data with model 3 cosmological pa- 
rameters. 

More precisely, taking into account the effective number 



Figure 4. The upper plot shows the WMAP best case (model 
1). The reconstructed matter power spectrum convolved with the 
SDSS window functions and evolved at redshift zero is given by 
the black diamonds connected by a piecewise continuous line. The 
small error bars reflect the quality of WMAP data and are con- 
structed with the square root of the diagonal elements of the full 
covariance matrix, which is used for the statistical analysis. The 
thin irregular curve is the mean reconstruction before convolution 
and the dash-joined triangles are the SDSS data with a maximis- 
ing bias. The lower plot is the SDSS best case (model 2) 



of degrees of freedom defined above, it is feasible to compare 
the probability of model 3 with the power law case, consid- 
ering both the WMAP and SDSS data. For a power law the 
probability is 3.9%, considering as free parameters h, Qt, 
O c , r, 6, n s and the amplitude of the fluctuations. While 
for model 3, with free parameters h, Qt, ^c, t, b and the 
effective contribution from the smoothed amplitudes that 
amounts to 45.2, the probability is 11.9%. 

In addition to a global goodness of fit, a statistical 
test to search for local features can be implemented. We 
combined the SDSS data and the SDSS-convolved WMAP- 
reconstructed power spectrum for model 3 to show more 
clearly where it is more likely to have deviations from a 
power law. The resulting power spectrum is depicted in 
Fig. [7] and its error covariance matrix and mean are given 



by 



£ _1 + 



f N 
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V 64 
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V 64 



b 2 



(6) 



(7) 



The feature at k ~ 0.05 /iMpc -1 stands out as a deviation 
from a power law at a significance of one sigma. Also note 
that the interesting eighth and ninth points are positively 
correlated at a level of approximately 20%. Similar detec- 
tions in other surveys suggest that this feature may be real 
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Figure 5. The conventions are the same as for the previous 
figure. In the lower plot, model 5, with a low baryonic content, 
presents oscillations completely out of phase with some of the 
SDSS features. In the middle, model 3, our best case, mimics 
particularly well the SDSS feature at k ~ 0.05 /iMpc -1 . Model 4 
at the top represents a case with slightly higher baryonic fraction 
and is also a good fit. 
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Figure 6. The mean (central curve) of our best case recon- 
structed matter power spectrum, model 3, is shown. The sur- 
rounding curves are obtained by summing and subtracting the 
square root of the diagonal elements of the covariance matrix. 
This is strictly legitimate only for an uncorrelated case, however 
here it is used to convey the magnitude of the uncertainties. The 
thick line is calculated from a primordial power spectrum given 
by k n "~ 1 , where n B sa 0.99. Most of the detected features do not 
seem to be the product of baryonic oscillations (for a possible ex- 
planation see e.g. Martin & Ringeval 2004). For convenience, the 
SDSS data are also depicted. 
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Figure 7. The power spectrum produced from a primordial 
power law with fc ns_1 , where n s ss 0.99 (smooth line) is depicted 
together with the combined spectrum from the SDSS data (data 
points with error bars) and the SDSS-convolved reconstruction 
for model 3 (see text). The SDSS feature at k ~ 0.05 /iMpc -1 
appears to be present at 1 sigma c.l.. 

and certainly encourage further searches in future experi- 
ments. 

The resolution of the released SDSS data does not al- 
low us to test the reality of various fine features that seem to 
be present in the reconstructed power spectrum. Along with 
progress in large scale structure measurements, a detailed 
understanding in forthcoming CMB experiments of possi- 
ble systematics in beam uncertainties and eventual non- 
gaussianities coming from foregrounds and point-sources 
residues would definitely increase the chances of detection 
of features in the primordial power spectrum. 



4 DISCUSSION 

In summary we have shown that a slight discrepancy in the 
baryonic fraction that arises from SDSS and WMAP data 
could be resolved if the primordial power spectrum is not 
fixed a priori in the form of a power law, but is constrained 
to satisfy CMB data when confronted with large-scale struc- 
ture information and letting key cosmological parameters 
vary. Our method is feasible due to a high resolution re- 
construction of the power spectrum from the CMB angular 
power spectrum, that is carried out for the first time in 
the full range considered, and provides evidence in favour of 
intrinsic features in the primordial power spectrum. In par- 
ticular a feature at k ~ 0.05 /iMpc -1 seems to occur in both 
the WMAP and SDSS data sets. This kind of technique of- 
fers an effective way to break the degeneracy between the 
baryon abundance and the shape of the power spectrum (in 
the conventional approach parametrized by the spectral in- 
dex, n s ) and relevantly brings also to a measure of the galaxy 
bias. 

A power spectrum with the observed features brings to 
a good fit of SDSS data without the necessity of a large 
baryon fraction. More generally, if the initial power law hy- 
pothesis is relaxed, deviations seem to induce a somewhat 
better concordance picture from the combination of CMB, 
large-scale structure, BBN and Hubble constant measure- 
ments. This could have important consequences for funda- 
mental physics. 
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